Astronomy <fc Astrophysics manuscript no. DzhaKuSta-R 


© ESO 2008 


September 2, 2008 





Wave instabilities in an anisotropic magnetized space plasma 

(Research Note) 

N.S. Dzhalilov 1 ' 2 ' 3 , V.D. Kuznetsov 2 and J. Staude 1 



00 

o 
o 

(N 
Or 
00 
(N 



6 



Astrophysikalisches Institut Potsdam (AIP), An der Sternwarte 16, D-14482 Potsdam, Germany, 

Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation of the Russian Academy of Sciences 

(IZMIRAN), Troitsk City, Moscow Region, 142190 Russia, 

Shamakhy Astrophysical Observatory of the Azerbaijan Academy of Sciences (ShAO), Baku Az-1000, Azerbaijan 
e-mail: JStaude@aip.de; Namig@izmirELii.ru; KVD@izmirEtn.ru 



Received ... ; accepted 1 July 2008 



ABSTRACT 



Aims. We study wave instability in an collisionless, rarefied hot plasma (e.g. solar wind or corona). We consider the 
anisotropy produced by the magnetic field, when the thermal gas pressures across and along the field become unequal. 
Methods. We apply the 16-moment transport equations (obtained from the Boltzmann-Vlasov kinetic equation) includ- 
ing the anisotropic thermal fluxes. The general dispersion relation for the incompressible wave modes is derived. 
Results. It is shown that a new, more complex wave spectrum with stable and unstable behavior is possible, in contrast 
to the classic fire-hose modes obtained in terms of the 13-moment integrated equations. 
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1. Introduction 

An almost collisionless, rarefied, hot, magnetized space 
plasma such as that of the solar corona is anisotropic and 
inho mogeneous, in part icular in the cross-field direction 
(see Aschwanden (2005) I. There have been observations of 
a thermal anisotropy of T±/Tu ~ 2 — 3 in the solar wind 
by both |Feldman et al. (1974)| and |Marsch et al. (1982)| of 
large heavy-ion thermal anisotropies (Tj_/T| > 100) by 
Kohl et al. (1998)land|Cranmer et al. (1999)1 of protons by 
Cranmer et al. (1999)| and of the coronal nole tempera- 



ture by Dodero et al. (1998) and Antonucci et al. (2000) 
An opposite ion temperature relati on T\\ > T± is also 



found in solar wind observations (see Marsch (2006)). Due 
to the anisotropy in the kinetic temperatures of protons 
and heavy ions, the corresponding partial pressures be- 
come anisotropic, in addition to the total thermal pres- 
sure such that p± =/= pii . It is now generally accepted that 
the observed large ion temperature anisotropies are related 
to the physical mechanism by which the solar corona and 
solar wind are heated (sec Hollweg & Isenberg (2002) and 
Marsch (2006)). 



In these circumstances, it is difficult to develop a tradi- 
tional hydrodynamical description of the plasma. We there- 
fore attempt to extend the MHD approximation by consid- 
ering the anisotropy of the magnetized plasma. We consider 
the large-scale wave peculiarities that can appear in a col- 
lisionless plasma. The large-scale plasma motions are usu- 
ally described by a fluid approximation, and the integrated 
moment equations derived from the Boltzmann-Vlasov ki- 
netic equation are used. If collisions between particles are 
rare and a strong magnetic field approximation is valid, 
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the usual MHD equations have to be replaced by other 
equations for which the fluid approximation is valid too. 
It was shown, e.g. by|Grad (1949)||Chew et al. (1956)[ and 
Rudakov & Sagdeev (1958) that for a collisionless plasma 
- mainly across a magnetic field - the fluid approach can be 
used. However, these so-called 13-moment equations can- 
not be used to describe a plasma of arbitrary anisotropic 
pressure. We therefore use the 16-moment equations. 



2. Basic equations and wave equations 

The 16-moment set of equations was used by 
many authors in different theoretical approaches, 
especially for modeling the solar wind (see 
IDemars fe Schunk (1979)] |Qlsen fc Leer (1999)[ p (1999) 



and Lie-Svendsen et al. (2001)p. The 16-moment set of 



transport equations for the collisionless plasma in the 
presence of gravity g b ut without magnetic di ffusivity is 
given as follows (see e.g. Oraevskii et al. (1985)): 



dp 
~dt 



pdivv = 0, 
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d_S± 

dB 

~dt + 



Pj| 

<i2 



P/ P P||-B 



(h-V)B 



Bdivv - (£? ■ V)t> = 0, divB = 0, 



(6) 
(7) 



where V = Vn + V_l,V|i = h(h ■ V), 



d 
dt 



dt 



(v ■ V), v = vn+ Vj_, h = 



B 

~B' 



(8) 



and S\\ and Sj_ are the heat fluxes along the magnetic field 
of parallel and perpendicular thermal motions. If the ther- 
mal fluxes are neglected, S± = and = 0, we obtain the 
equations describing the laws of the change in longitudinal 
and transverse thermal energy along the trajectories of the 
plasma (the left-hand parts of Eqs. ([3]) and These so- 
called "double-adiabatic" invariants and Eqs. (|T|), and 
also form a closed system of eq uations, the CGL (C hew- 
Goldberger-Low) equations (see Chew et al. (1956) ). By 
using the CGL-equations, we would however obtain incom- 
plete equations instead of Eqs. (0 [5]). This is because, by 
deriving the CGL equations, authors so far omitted with- 
out proof the third moments of the di stribution function 
and therefore the thermal fluxes (se e Chew et al. (1956) 
and Baranov & Krasnobayev (1977) ). The equations de- 
rived for the 16-moment set, in our case Eqs. HH7]), in- 
clude the thermal fluxes; they are more complete, and the 
CGL equations cannot be derived from these equations 
as a special case. For simplicity, we assume that the ba- 
sic initial equilibrium state of the plasma is homogeneous 
(g = 0, and the quantities vq, Po,P±o,P\\o, Bq, 5j_o, and S\\q 
are constant). Equations {T]) to automatically satisfy 
this equilibrium state with non-zero thermal fluxes. We con- 
sider small linear perturbations of all physical variables, 
for example pressure in the form p = po + p'(r, t), where 
p'(r,t) ~ expj(fc • r — tot), uj is the wave frequency, and 
k is the wave number. For the perturbations, we obtain 
linear wave equations. Even if we insert zero initial heat 
fluxes = Sj_o = 0, the perturbations of these functions 
will never become zero: S'^ ^ 0, S' ± ^ 0. Using the 16- 
moment equations, we should derive more reliable results 
about the wave properties in an anisotropic plasma than 
with the CGL equations based on the 13-moment equa- 
tions. 

In the presence of an external magnetic field, the initial 
collisionless heat fluxes should be defined by solutions of the 
kinetic equations. We should, however, use some appropri- 
ate estimate as a parameter. The heat flux functions should 
be estimated by taking the thermal energy density of the 
electrons multiplied by the particle stream speed along the 



magnetic field u : S 



|o 



k B T\\u S = j8uoP\\. Hollweg 



(1974, 1976) provided some estimates of the correction pa- 
rameter 6 (a in these papers) by assuming realistic shapes 
of electron distribution functions and comparing the results 
with space observations. We note that 8 depends on the 
magnetic field. In the range of B = 0.1 — 100 G, the esti- 
mates give 5 « 4 — 0.1. In the same way, S±o ~ j5u p±, 
and we define the parameter 7 = (3/4)(5u / c || ■ We note 



that Marsch & Richter (1987) quoted values of 7 measured 
in the solar wind. 



We introduce dimensionless parameters and note that 
the indices " 0" of physical parameters are omitted for sim- 
plicity: 



a = 



P± 
P\\ ' 



P\\ 



Sn = 



,0 



Si 



B 2 
47rp|| 



P||C|| 



P±c\\ 



cos 



(9) 



(10) 



where (f> is the angle between wave vector and magnetic 
field, and the indices || and _L correspond to the values 
of the parameters along and across the magnetic field, re- 
spectively. We note that our (3 is defined to be inversely 
proportional to the more commonly used plasma beta, 
(3 = 2//3 p i asma . With the parameter 7 defines above, we 



have 5*11 = S 



7- 



In analogy with the usual MHD equations used e.g. 
by Somov et al. (2007)] there are two independent wave 
branches in the plasma: waves that do not compress the 
plasma (div v = 0) and waves that compress the plasma 
(divw 7^ 0). We restrict ourselves to the incompressible 
wave modes. 

After inserting into the wave equations the condition of 
incompressibility (k ■ v) — and p' = 0, we obtain as usual 
the parametric dispersion equation. This is a polynomial 
equation of 6 th order in the frequency of oscillations. For 
the parameter Z = w/(c||fc||), this equation can be written 
in the form 

c 6 Z 6 + c 5 Z 5 + c 4 Z A + c 3 Z 3 + c 2 Z 2 + Cl Z + c = 0, (11) 
c - 3 (1 - a) [a 2 (l - 1) + a(l + I) + 2 {fi - I)] , 
ci = 2 7 [/3(a-2) + 3a(l-2Z) + a 2 (4;-3) + 2Z] , 
c 2 = 4Z + 2/?(4a- 5) + 2a 2 (l + 3Z) + a 3 (Z- 1)- 
-a(l + 110, c 3 = -2 7 [j3(a - 2) + a 2 (2 1 - 1) - 2 a] , 
c 4 = 6 1 + a(l - 3) - 2 la 2 + 2 /3(2 - a), 



C5 



27(2ia - 21 — a) , c 6 = a + la-4l. 



All coefficients are real and, consequently, all solutions are 
real or conjugate complex. In the usual isotropic MHD case, 
only Alfven waves with uj 2 = k 2 v\ are present, the phase 
velocities of which are equal to each other in both directions 
with respect to the magnetic field. Instead of assuming that 
Z 2 = P in the isotropic MHD, we determined now the 6 th 
order Eq. (fTTj) in the anisotropic case. With the heat fluxes 
for which 7^0, odd nonzero coefficients 01,03, and C5 
generate wave propagation velocities that depend on the 
direction of the magnetic field. 



3. Limiting and special cases 

In this section we investigate the most important cases of 
Eq. (TTTj) which can be solved analytically. 
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3.1. Parallel propagation 

In the case I = 1 or k = fc||, the six roots of Eq. pi[) arc 
simple. The first pair of roots corresponds to a pair of stable 
modes 



(12) 



The properties of these modes do not depend on the mag- 
netic field. They are isotropic in terms of the direction of 
magnetic field. This implies that stable waves propagate 
along and against the magnetic field with the same phase 
velocity. In these modes, p' = 0, v' = 0, B' = 0, pj, = 0, 

and 5*11 =0. However, p' ± ~ S' ± ^ 0. The restoring force for 
these modes is therefore Wp± . Motions of particles across 
the magnetic field cause a perturbation of the plasma pres- 
sure p_L, but this is compensated by the generation of a heat 
flux. At the same time no hydrodynamic motion is produced 
by these modes, v — as pj, — 0. We therefore have unusual 

thermal waves: they are stable waves (Im(w) = 0) that 
propagate (Re(w) 7^ 0) with the parallel sound speed cn. 
Analogous to sound waves, they are not dispersive modes. 
The modes of this branch are named "isotropic thermal" 
waves. 

The second pair of roots corresponds to 



w/fc||C)| = ±\Ja + (3-1, 



(13) 



which are conventional isotropic fire-hose waves. The waves 
become unstable if a + (3 < 1 or p» > p± + 2p mag . The fire- 
hose instability therefore disappears for (3 > 1 or a > 1. 
The dispersion relation (fT3|) passes into that for the usual 



(14) 



isotropic Alfven waves lo 2 = k 2 v\ if a = 1. 
The third pair of solutions corresponds to 



2— = ^±^ + 12^. 



The instability condition (a— l)/(2 — a) > 7 2 /12 is obeyed 
if 1 < a < 2. For 7 / 0, the unstable modes begin to 
propagate, Re(cj) 7^ 0. For 7 = 0, stable waves outside the 
region 1 < a < 2 travel along and opposite to the direction 
of the magnetic field with the same phase velocity. With 
7^0, the stable waves for which (a — l)/(2 — a) < 7 2 /12 
become anisotropic. Retrograde waves travel more rapidly 
than prograde waves. These branches are named thermally 
"anisotropic" waves. With increasing 7, the prograde waves 
propagate more slowly, but retrograde modes travel faster. 
There is a crossing of these branches on the axes of a for 
7 = at V p h = 0, and for 7 > at V p h < 0. Between the 
crossing points of the branches and a — 2, an instability 
arises. In contrast to the fire-hose and the thermal modes, 
these modes in the instability region are traveling, Re(cj) 7^ 
0. They do not depend on (3. For these waves 



n 

S'n 



2(1 



B' 



EL 
Pi 



1 - a 2 B> 



\-2otB 



1-a B' 
l~2a~B' 



These modes are generated by pressure anisotropy: if a = 1, 
they disappear. 



3.2. Oblique propagation 
3.2.1. Strong magnetic field 

The limit f3 ^> 1 produces solutions of Eq. lfTTj) that 
are similar to those in the parallel propagation case. In 
this case, the fire-hose waves become stable with high 
phase velocities, Z 2 ~ 0{(3). Other isotropic thermal and 
thermally anisotropic waves remain unchanged. 



3.2.2. Weak magnetic field 

In the limit (3 <C 1, the first pair of solutions of Eq. (Ti"Tj) sim- 
t — 1, which is the same solution as that for 



plifies to Z w 

the parallel fire-hose waves. The other thermal modes are 
described by a polynomial equation of 4 th order. However, 
for the limit 7 — > 0, we obtain 



Z 2 sa (-b*±VD)/(2a*), D 



4a*c* 



(15) 



where c* = 3(1 — a)[a(l — I) — 21], a* = a + l(a — 4), 
and 6* = a 2 (l — I) — 4a(l + I) + 101. These solutions can 
easily be investigated. In the region < a < 1 for a closed 
area (a, I), we have D < and, therefore, Re(Z) ^ and 
Ixa(Z) 7^ 0. This implies that both thermal waves become 
unstable, and D > if a > 1. In this case for unstable 
modes, Re(Z) = and instability disappears if a > 2. The 
thermal wave branches have crossing points, where their 
phase velocities coincide. For example, this occurs if / = 1 
and a = 0.5. 



3.2.3. Isotropic propagation 

We consider the important special case of arbitrary [3 
and 7 = 0. Even though the absence of fluxes, S± = 
and S'n = 0, is far from reality, this simplified 
case was investigated u sing the 13-moment equations 



(e.g. Kato et al. (1966) 



Baranov fc Krasnobayev (1977) 



and Kuznetsov &: Oraevskii (1992)[ ). Substituting 7 = 6 



into Eg. (jlip. we obtain a cubic equation for Q = Z 
c 6 C 3 + c 4 C 2 + c 2 C + c = 0. 



(16) 



We have symmetrically only three pairs of solutions, Z = 
±\/C- The analytical solutions C11C2, and £3 allow us to 
investigate in significant detail the dependence of the so- 
lutions on the parameters a, f3, and <fi. For a small devia- 
tion of the propagation angle (</> ^ 0), the situation dif- 
fers significantly from that of the parallel propagation. All 
three modes interact, and this interaction occurs in two 
ranges of a: a > a c and a < a c . The critical value of 
a = a c = + is the singular point for cq = 0. 

With increasing propagation angle, the interaction domains 
expand. In these domains, the growing rates also become 
larger and we obtain a mixture of modes — a turbulent 
wave motion. 

3.2.4. Anisotropic propagation 

To investigate the role of the thermal parameter 7, we 
solved the 6 th order polynomial equation, choosing more 
realistic values of 7 < 1 . In Fig. ([1} , the phase velocities and 
instability rates are shown. We observe that 
a) waves with positive and negative velocities arc different 
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Fig. 1. Phase velocities V p h — Re(Z) of the three modes as a function of a for fixed parameters = 1, <fi = 7r/4, and 
7 = 0.1 (left picture); instability rates Im(Z) of the three modes (right picture). The labels at the curves correspond to 
modified fire-hose waves (/), to isotropic thermal waves (i), and to thermally anisotropic waves (a). 



and all three wave branches become coupled; 

b) maxima of the instability rate strongly depend on 7. 

For some parameters, all three modes become unstable. 

3.3. Classic fire-hose modes 

To compare our results with the classic fire-hose modes 
based on the CGL equations, we use the initial set of equa- 
tions, omit Eqs. ([5]) and ©, and substitute into the others 
Sn = S± = 0. These CGL equations for incompressible 
waves provide two pairs of solutions: 

z = ±iVA, A ^(l^Mlp_Jll±W^l. (17) 

a(l + l)-4l 

The fire-hose instability condition is that A > 0. In the 
parallel propagation case (I = 1), this condition is the most 
familiar case of a + (3 < 1. We obtain the same result for 
I 7^ 1, if f3 = 0. The instability condition is more complete 
for the oblique propagation case if [3 > 0, but this strongly 
differs from our results based on the 16-moment equations. 
For the classic fire-house instability, we always have that 
Re(w) = 0. In our case, Re(w) 7^ in most cases due to 
the coupling of these modes with other thermal modes. For 
the classic modes in the relations, for example, = 
niB'/B and p' ± /p± = ri2B'/B, the coefficients are n\ = 
—2, ri2 = 1. In our case, these coefficients are complete 
functions of all parameters. The main difference in our case 
compared to that for the CGL equations is the appearance 
of two additional thermal branches, even if 7 = 0. 



4. Conclusion 

To investigate the peculiarities of large-scale wave motions 
in a collisionless magnetized plasma, we have applied the 
16-moment transport equations, derived as integrated mo- 
ments of the kinetic equations. In earlier similar attempts, 



the 13-moment equations were used. However, these equa- 
tions exclude without any reason the thermal fluxes and 
are therefore incomplete. 

Anisotropy is the main feature of a collisionless plasma 
with a strong magnetic field. In the present study, the pres- 
sure anisotropy was described by the parameter a and the 
heat fluxes by 7. By assuming that 7 = we were unable to 
derive the 13-moment equations, or by assuming that both 
a = 1 and 7 = we did not obtain the isotropic MHD 
case. The 16-moment equations were in principle differ- 
ent equations. Using these equations, we illustrated that a 
wide unstable and stable wave spectrum in the collisionless 
anisotropic plasma was possible, even in the incompress- 
ible approximation. If 7 =/= (heat fluxes are present), the 
waves propagated along and against the magnetic field at 
different speeds. This behavior differed from that of the 
usual isotropic MHD case. The coupled wave spectrum, in- 
cluding modified fire-hose modes, strongly depends on the 
magnetic field value (parameter /?), pressure anisotropy pa- 
rameter a, heat flux parameter 7, and wave propagation 
angle 4> with respect to the magnetic field. The deduced in- 
stability increments are rather large. We have derived the 
general instability condition for incompressible waves. 
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